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ABSTRACT 

Recently, eigenmodes of rotating fluids, namely inertial modes, have received much attention 
in relation to their destabilization when coupled to gravitational radiation within neutron stars. 
However, these modes have been known for a long time by fluid dynamicists. We give a short 
account of their history and review our present understanding of their properties. Considering 
the case of a spherical container, we then give the exact solution of the boundary (Ekman) layer 
flow associated with inertial r-modes and show that previous estimations all underestimated 
the dissipation by these layers. We also show that the presence of an inner core has little 
influence on this dissipation. As a conclusion, we compute the window of instability in the 
temperature/rotation plane for a crusted neutron star when it is modeled by an incompressible 
fluid. 

Subject headings: Hydrodynamics; rotating stars; neutron stars 



1. Introduction 

Recently much work has been devoted to the 
study of the rotational instability of neutron stars 
resulting from a coupling between gravitational 
radiation and the so-called "r-modes" of a rotat- 
ing star Andersson (1998); Friedman and Morsink 
(1998); Lindblom et al. (1998); Kokkotas and Ster- 
gioulas (1999). Such an instability may indeed 
play a key role in the distribution of rotation pe- 
riods of neutron stars as well as it may be an im- 
portant source of gravitational radiation. 

In this paper, we shall first clarify a point of 
history concerning "r-modes" which are in fact a 
special class of inertial modes; we shall then review 
their singular properties which have been clari- 
fied only very recently in Rieutord and Valdettaro 
(1997), Rieutord et al. (2000a) and Rieutord et al. 
(2000b). The last section will present the ana- 
lytical derivation of the damping rate of inertial r- 
modes in a neutron star with a crust and/or a core 
through the boundary layer analysis within the 



framework of newtonian theory. We conclude on 
the stability of crusted neutron stars when mod- 
eled by an incompressible viscous fluid in a rotat- 
ing sphere. 

2. A short point of history 

The very first work on rotating fluids oscilla- 
tions which are presently known as inertial modes 
dates back to Thomson 1 (1880) who analysed the 
case of a fluid contained in a cylinder. How- 
ever, another impetus to the study of these oscil- 
lations was given soon after by Poincare's (1885) 
work on the stability of rotating self-gravitating 
masses, a work applied to MacLaurin spheroids 
by Bryan (1889) 2 and later continued by Car- 
tan (1922) who christened the equation of iner- 
tial modes as "Poincare equation" . In these stud- 
ies, however, the effect of rotation is combined to 
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2 but see the recent rederivation by Lindblom and Ipser 
(1999). 



the one of gravity through (for an incompressible 
fluid) surface gravity waves. In fact, except for 
the work of Thomson, investigations on the os- 
cillations specific of rotating fluids seem to have 
started with the work of Bjcrknes et al. (1933) 
where they are called "elastoid-inertial oscilla- 
tions" since conservation of angular momentum 
makes axis-centered rings of fluid behave elasti- 
cally; but see Fultz (1959) or Aldridge (1967) for 
an account on this part of history. In the sixties, 
much work has been devoted to these oscillations, 
mainly by Greenspan who introduced the termi- 
nology of "inertial oscillations" . The presently 
used denomination "inertial modes" has been "of- 
ficially" given by Greenspan's book (Greenspan 
1969). 

However, inertial modes are somewhat too gen- 
eral for applications in some specific domains like 
atmospheric sciences. In this field indeed motions 
are essentially two dimensional and inertial modes 
may be simplified into the well-known Rossby (or 
planetary) waves. 

The introduction of r-modes by Papaloizou and 
Pringle (1978) was quite unfortunate since they 
associated eigenmodes of rotating fluids with a 
very special class of inertial modes, namely purely 
toroidal inertial modes. This lead following au- 
thors to introduce weird names such as "hybrid 
modes" or "generalized r-modes" (Lockitch and 
Friedman 1999) for describing the general class of 
inertial modes. We therefore encourage authors 
to use, as fluid dynamicists, inertial modes unless 
they discuss the very specific r-modes. 

3. The present theory of inertial modes 

Inertial modes are a class of modes of oscillation 
of rotating fluids which owe their existence to the 
Coriolis force. This force of inertia has indeed a 
restoring action on perturbations of rotating fluids 
since it insures the global conservation of angular 
momentum. These modes have many properties 
similar to those of gravity modes of stably strati- 
fied fluids (Rieutord and Noui 1999). 

The dynamics of inertial modes may be ap- 
preciated when all other effects are suppressed: 
no compressibility, no magnetic fields, no gravity, 
etc... only an incompressible inviscid rotating (like 
a solid body) fluid. In this case, perturbations of 
velocity 5v and pressure SP obey 



dt 



+ 2n x 5v = -V5P, V • <5v = (1) 



where ft is the angular velocity of the fluid. Con- 
centrating on time-periodic oscillations and choos- 
ing (2fl)~ 1 as the time scale, (1) can be written 



iuju + e z x u 



-Vp, 



V-u = 



(2) 



with non-dimensional variables; uj is the non- 
dimensional (real) frequency. When the velocity 
u is eliminated in favor of the pressure perturba- 
tion p, one is left with 



Ap- — - 

uu 2 dz 2 



(3) 



which is known as Poincare equation since Car- 
tan (1922). This equation is remarkable in the 
fact that it is hyperbolic spatially since \ui\ < 1 
(Greenspan 1969). As the solution of (3) must 
meet boundary conditions, namely u • n = 0, we 
see that inertial modes are solutions of an ill-posed 
boundary value problem 3 . This property means 
that, in general, inertial modes are singular; in 
other words they cannot exist physically if the 
fluid is strictly inviscid. These properties are de- 
tailed in Rieutord and Valdettaro (1997) and Rieu- 
tord et al. (2000a); to make a long story short, one 
may summarize the situation as follows. 

Let us first recall that in hyperbolic systems, 
energy propagates along the characteristics of the 
equation. For the Poincare problem, these are 
straight lines in a meridional plane. A way to ap- 
proach the solutions of this difficult problem is to 
examine the propagation of characteristics as they 
reflect on the boundaries. They define trajecto- 
ries which depend strongly on the container. Let 
us therefore concentrate on the case of a spherical 
shell as a container; this configuration is relevant 
for neutron stars with a central core due to some 
phase transition of the nuclear matter (see Hacnsel 
1996). In this case, it may be shown that charac- 
teristic trajectories generically converge towards 



3 Such boundary conditions eliminate any distorsion of the 
surface due to the fluid motion; for a free surface, these 
distorsion are surface gravity waves (see Rieutord 1997) 
but their inclusion (in order to be more realistic) would 
not modify the ill-posed nature of the problem. 
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Fig. 1. — The kinetic energy distribution in a 
meridional plane of an inertial mode in a spher- 
ical shell associated with an equatorial attractor. 
A co-existing polar attractor is also slightly ex- 
cited. The mode is axisymmetric with equato- 
rial symmetry. Stress-free boundary conditions 
have been used on both shells; this solution was 
computed with an Ekman number of 2 10 -9 and 
required 1300 spherical harmonics and 450 ra- 
dial grid points (Gauss-Lobatto). The ratio of 
the inner radius to the outer radius is r\ = 0.35. 
lu = 0.2429 and the damping rater = — 6.26xl0~ 4 
are given in dimensionless units as equation (2). 

attractors which are periodic orbits. It may be 
shown (Rieutord et al. 2000a) that in this case, the 
associated solutions are singular, namely the ve- 
locity field is not square- integrable. However, iner- 
tial r- modes are still solutions of the problem since 
they meet the boundary conditions (u r = 0); in 
fact they are the only regular (square-integrable) 
solutions of the Poincare problem in a spherical 
shell. In a more mathematical way, we may say 
that the spectrum of eigenvalues of the Poincare 
problem in a spherical shell is empty except for the 
inertial r-modes. In this sense, these modes are 
quite exceptional. This situation occurs because 
there exists no system of coordinate in which the 
dependent variables of the Poincare equation can 
be separated. This is a consequence of the con- 



flict between the symmetry of the Coriolis force 
(cylindrical) and the geometry of the boundaries. 
Thus, when this constraint is relaxed, like in the 
case of a cylindrical container, regular solutions 
exist and a dense spectrum of eigenvalues appears 
in the allowed frequency range, namely [0, 2f2]. In 
the case the container is a full sphere, attractors 
also disappear and eigenmodes exist; they are also 
related to a dense spectrum of eigenfrequencies. 
In this case, Poincare equation is exactly solvable 
(Greenspan 1969). 

However, real fluids have viscosity (y) and 
equation (2) should be transformed into 

Au + e z x u = -Vp + EAu, V • u = (4) 

where A is the complex eigenvalue and E = 
v/2Q,R 2 is the Ekman number (R is the outer 
radius of the shell) . 

Using no-slip (u = 0) or stress-free boundary 
conditions, (4) yields a well-posed problem. Yet, 
the singularities of the associated inviscid solu- 
tions show up through the existence of shear lay- 
ers. As shown by fig. 1, the shape of inertial modes 
is deeply influenced by the underlying singularity 
of the inviscid solution. We have shown (Rieutord 
et al. 2000a) that these shear layers are in fact 
nested layers with different scales for their inner 
part scales as i? 1 / 3 and their outer part seems to 
scale with E 1 / 4 . Because of these internal shear 
layers, these modes are strongly damped. 

We therefore see that according to whether a 
neutron star has a central core or not, the damp- 
ing of inertial modes will be extremely different. If 
there is a central core, the only regular modes are 
the inertial r-modes which will be by far the least 
damped; if there isn't any core then a dense spec- 
trum exists (Greenspan 1969; Lockitch and Fried- 
man 1999) but inertial r-modes remain the most 
unstable because of their simple structure. 

4. Damping of toroidal inertial modes in a 
sphere or a shell 

We shall now give the expression of the vis- 
cous damping of inertial r-modes when one of the 
boundary is solid therefore when the dissipation 
is due to Ekman boundary layers; we shall thus 
complete the works of Bildsten and Ushomirsky 
(2000) and Andersson et al. (2000) by giving the 
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rigorous estimate of the damping rate; the method 
which we use here has been outlined in Greenspan 
(1969). 

The damping rate is given by: 

(•) 

where (V : u) 2 stands for the squared rate-of- 
strain tensor (see below). The velocity field of 
r- modes is 

ug = Ar m (sm6) m ^ 1 sin(m0 + u m t) 
u v = Ar m (sin 9) m ^ 1 cos 9 cos(m(j) + uj m t) 

The kinetic energy integral may be evaluated 
explicitly 



Now we need the expression of the square of the 
rate-of-strain tensor Sij = diVj + djVi in spherical 
coordinates, viz 



(V : u) 2 = s 2 rr + s 2 e + s 2 + 2(s 2 r0 + s 2 + 



Since the radial derivatives dominate, this ex- 
pression reduces to the contribution of the tan- 
gential stresses. Using the scaled coordinate, ( = 
\r — rb\/VE, we have 



(V:u) 2 = ^ 



/ dug 



+ 



We now set p = cos 9 — ui and q = cos 9 + u). We 
thus get 



/ 



u 2 dV = irA 2 (1 - r, 2m+3 ) + !m' 
v ' ; m(2m + 3)!! 



(6) 



where r\ is the ratio of the radius of the inner 
boundary to the one of the outer boundary. 

The dissipation integral need more work if one 
of the boundary is no-slip. In this case, dissipation 
is essentially coming from the Ekman layers and 
thus we need to derive the flow in these layers. The 
method has been given by Greenspan from whom 
we know that the boundary layer correction u, is 
related to the interior solution u by 



uq + iu v = — (ug + iu v ) r = rb e 



(7) 



where (, is the radial scaled variable (r — r^/y'E 
with rfc as the radius of the boundary (1 or rj). 
The complete solution is then u = u + u; setting 
(3 = cut + nrup we have 



Ug + iu v = 



Ar m (sin 1 



\m-l 



2i 



x {(l-cos60e l/3 -(l + cos6>)e- 4/3 } 



from which it follows that 



(V : u) 2 = A2r b m ^) 2m 2 j (1 + cos^ble-f^ 

+ (l-cosfl) 2 |< ? |e- c \ /2 i^ 
-2 sin 2 6 3?e (e 2 ^C(V*9+^P)^ J 

Integrating over the (/^-variable yields 



2tt 



(V : ufdip = 



it A 2 r 2m (sin 6) 
E~ 



2m-2 



x {(l + cos0) 2 |p|e- c V 2 l^+(l-cos0) 2 | g |e- c V 2 l^} (8) 
We now integrate over the radial variable: 



T) Jo 



(V : ufdip r 2 dr 



7rA 2 (sin6») 



2m -2 



E 



f 



,2m+2 



f«)dr 



with /(C) = (1 + cos9) 2 \p\e-t^\ + (1 - 
cos#) 2 |(7|e~''V /2 ~i<?l; since r = 77 + VE( or r = 
1 — ^/EC > according to which side of the integral is 
chosen, it turns out that 



Ug + lUu 



Ar™(s\n9)' 
2i 



-ij3—C,^/i cos d — iuj 



- {(1 + cos <9)e 

- (1 - cos6')e^ CVicos0+ ^| 



(V : u) 2 dif r 2 dr 



T) JO 



r-OO 

K / (l + cos6») 2 |p|« 
Jo L 



/2\p\ 



4 



Outer B.C. no-slip 



stress- free 



Inner B.C. 

no-slip 
stress-free 



l + ?7 



2m+2 



2m+2 



1 



T] 

No Ekman layer 



Table 1: Expression of the P(rj) function according 
to boundary conditions. 



+ (1 -cos6) 2 \q\e- c V2M 
K 

71 



(l + cos0)Vl cos ^-^l 



+ (1 - cos6») 2 v /|cos6» + w| 



with K = 7 r J 4 2 (sin6') 2m - 2 P(r7)/VS, where P{rj) is 
a function depending on the boundary conditions 
(see table 1). Finally integrating over 9, we find 



/ 



(V : uf dV = 2 ^!£.<">I.. 



>2E 



(9) 



with 



I m = / (l + cosfl)V|cosfl- 
Jo 



- uj \ sin 



2m- 1 



(10) 

Finally grouping (6) and (9), we find the damp- 
ing rate 



7 = 



m(2m + 3)!! 
"2 m + 3 / 2 (m + l)! 



Q( V )l m VE (11) 



where Q(r?) = P(rj)/ (l - 77 2m + 3 ). 

For the cases m = 1 and m = 2 we evaluated 
the expression of l m , viz 



Xi = 



V2 
35 



(3 5 / 2 + 19) 



2 \ 11/2 3401 + 2176^ 



5 x 7 x 9 x 11 

Other values are computed numerically and given 
in table 2. 

The values given by (11) may be compared to 
other derivations, in particular that of Greenspan 
(1969) for m = 1 who finds -y/y/E = -2.62/V^ = 
— 1.8526 ! For m = 2, a direct numerical calcu- 
lation, similar to that of Ricutord and Valdcttaro 



m = 1 



m 



m 



m = 4 



J m 1.3976 0.80411 0.58075 0.46155 
^ -1.8526 -2.4876 -3.0318 -3.5339 



Table 2: Values of the first integrals T m and the 
corresponding values of the scaled damping rates; 
this latter value equals that of Greenspan (1969) 
when multiplied by \J~2~ because of our choice of 
the time scale. 

(1997), gives -2.482VS at E = 10~ 8 which is in 
good agreement with the analytical formula. 

5. Application to neutron stars and con- 
clusions 

Let us now apply these results to the case 
of rapidly rotating neutron stars. We take the 
viscosity from Bildsten and Ushomirsky (2000), 
v = 1.8 //T| m 2 /s where / is a dimensionlcss 
parameter taking into account the different trans- 
port mechanisms in the fluid (supcrfluid phases 
for instance) and Ts is the temperature in 10 8 K 
unit. Using a radius of 12.53 km and an angular 
frequency of 27rxlkHz, we find an Ekman num- 
ber ~ 10~ 12 which is indeed very small and thus 
boundary layer theory applies. 

We may now estimate the charateristic time 
scale for the damping of the m = 2-mode. We 
find 



26.7s 



(- 



R 



yff \ 10 km/ V v. 



1 kHz 



1/2 



(12) 



which is a somewhat smaller value than the pre- 
vious estimate of Bildsten and Ushomirsky (2000) 
and Andersson et al. (2000) who find a charac- 
teristic time of 100 s and 200 s respectively. Our 
disagreement with these authors comes from their 
approximate evaluation of the boundary layer dis- 
sipation and from the resulting functional depen- 
dence with respect to mass and density. Let us 
first evaluate the damping rate according to Lan- 
dau and Lifchitz (1989); it turns out that 
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Fig. 2. — Curves of critical angular velocity, nor- 
malized by VLk = | \/7rG'p, for different models. 
The solid line shows the result of the present work, 
the dashed-dottcd one is that of Andersson et al. 
(2000) and the dashed one is for Bildsten and 
Ushomirsky (2000). The dotted line is the criti- 
cal curve for a non-crusted star. No core has been 
included (r] = 0). 



27 = -(- 



1/2 



(J^u 2 smffffldip) (r = 1) 



(13) 

where we used our non-dimensional units. Since 
the radial dependence of the modes is in r m and 
co = l/(ra + 1), we easily find that 



7 



2m + 3 
2y/2m + 2 



\[~E 



When this expression is applied to the m = 2- 
mode, we find that 7 = — 1.429a/B which is a fac- 
tor 1.74 weaker than the correct result. 

If we use, as previous authors, a step function 
for describing the density difference between that 
of the crust and the mean density, we find that the 
damping rate reads 



lEk = -2.4876^^ 212 = -0.0346^ ^PV 1 
P P T 8 

(14) 

where pb is the density of the fluid just below the 
crust and 12* = Cl/^/wUp. 

Our calculation therefore shows that the win- 
dow of instability in the 12, T plane is smaller than 
previously estimated for crusted neutron stars. 

Considering a 1.4 M neutron star with a ra- 
dius of 12.53 km as a test case, the growth rate of 
the mode due to gravitational radiation is "f gw = 
0.658s _1 12* (we use the expression given in Lind- 
blom et al. 1998); although, it is not relevant for 
an incompressible fluid, we take into account the 
damping rate due to bulk viscosity in order to ease 
comparison with previous work; from Lindblom 
et al. (1999), we find lbuik = -2.2 10~ 12 s" 1 T 9 6 12 2 . 
From (14), we have 7^ = -1.5310- 3 s" 1 nl /2 /T g 
where we took p b = 1.5 10 17 kg/m 3 ; solving the 
equation 

Igw + lEk + iBulk = 

for different values of the temperature yields the 
curves displayed in figure 2. 

As expected, we see that the window of instabil- 
ity narrows compared to Andersson et al. (2000): 
for a given temperature, the critical angular veloc- 
ity raises by <~10% typically. 

Another interesting conclusion of this work is 
that the presence of a solid inner core does not 
change the damping rates very much unless its 
radius is close to unity. The reason for that is 
to be found in the shape of the inertial r-modes 
whose amplitudes are concentrated near the outer 
boundary. Therefore, the rotating instability of 
rapidly rotating stars is quite insensitive to the 
presence of a solid core and more generally to any 
phase transition which does not occur close to the 
surface. 
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goulhon for drawing my attention on these ques- 
tions and for helpful discussions. I am also very 
grateful to Ian Jones for his note about the density 
profile used in models of neutron stars in previous 
work. Part of the calculations have been carried 
out on the Nec SX5 of IDPJS at Orsay and on the 
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gratefully acknowledged. 
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